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Abstract 

We consider a hyperbolic conservation law posed on an (N + l)-dimensional spacetime, 
whose flux is a field of differential forms of degree N. Generalizing the classical Kuznetsov's 
method, we derive an L 1 error estimate which applies to a large class of approximate solutions. 
In particular, we apply our main theorem and deal with two entropy solutions associated with 
distinct flux fields, as well as with an entropy solution and an approximate solution. Our 
framework encompasses, for instance, equations posed on a globally hyperbolic Lorentzian 
manifold. 

1 Introduction 

This paper provides a general framework leading to error estimates for hyperbolic conservation 
laws posed on an (N + l)-dimensional manifold M, referred to as a spacetime and, in particular, 
leading to a sharp estimate for the difference, measured in the L 1 norm, between an exact solution 
and an approximate solution. The present paper can be regarded as a generalization to manifolds 
of a contribution by Bouchut and Perthame [4], who recast in a concise form the pioneering works 
of Kruzkov and Kuznetsov (6HZH8) for hyperbolic conservation laws posed on the flat (Euclidian) 
spacetime. We are thus interested here in extending these results to conservation laws defined 
on manifolds, and develop a physically more realistic setting when geometrical effects are now 
taken into account. 

Motivated by the case of the shallow water equations on the sphere, the theory of hyperbolic 
conservation laws on manifolds has been developed in recent years by LeFloch together with col- 
laborators. In particular, well-posedness results have been obtained in l3ll9l [Til , and convergence 
results for finite volume schemes in (TJ [2), while an error estimate in the case of a Riemannian 
manifold was derived in flOl . For further results on the well-posedness theory, we also refer the 
reader to contributions by Panov in [12 . 13 1 and, on the finite volume schemes, to the earlier work 
0. 

Recently, in ITU , LeFloch and Okutmustur introduced a framework based on differential 
forms and dealt with conservation laws defined on an (N + l)-dimensional manifold M. In their 
formulation, the flux of the equation is given by a field of N-f orms, rather than by a vector field 
as was the case in earlier works. The formulation based on N-f orms is geometrically natural in 
that only minimal assumptions on the geometrical structure are required in order to establish the 
well-posedness of the initial value problem. 
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Considering hyperbolic conservation laws posed on a spacetime, we derive here a coordinate- 
free error estimate in the L -norm, which involves, on one hand, an exact weak solution satisfying 
all entropy inequalities and, on the other hand, an approximate solution for which the entropy 
dissipation measures contain some (positive) "error terms". Our assumptions on the manifold are 
minimal and, in particular, no metric is needed and no volume form is a priori prescribed. Indeed, 
our method only requires a suitable "hyperbolicity condition", and the prescription of a certain 
family of mollifiers (similar to mollifiers that may be determined from the distance function when 
a Riemann metric is given on the manifold). Still, we observe below that, to the flux-field defining 
the conservation law, one can associate a natural choice of volume form defined on the manifold. 
In addition, from our main theorem we deduce, on one hand, the L 1 -contraction property of the 
semigroup of entropy solutions and, on the other hand, a sharp estimate comparing together the 
solutions of two conservation laws with different flux fields. 

In the second part of this paper, we investigate the case of a conservation law defined on 
a Lorentzian manifold. Here, the presence of a metric allows us to refine our estimates. For 
instance, a suitable family of mollifiers is provided by the metric structure of the manifold, and 
the error terms can be given in a more explicit way, since they are naturally written using the 
metric. More importantly, the introduction of a metric allows us to consider second-order error 
terms in the approximate solutions. This allows us to apply our main result further, and derive 
error estimates for a nonlinear diffusion model. 

We may summarize the main difficulties overcome in this paper as follows. First of all, since 
the conservation law under consideration is posed on a non-flat manifold, the geometry of that 
manifold must be taken into account in, for instance, the formulation of approximate schemes 
and the analysis of their convergence. It has been pointed out that geometric effects occur which 
change the qualitative properties of solutions [1]. Spacetimes and, in particular, Lorentzian 
manifolds may not be "invariant by translation" in the time direction, and the time and spatial 
geometries are intertwined, giving rise to phenomena not present in the Euclidian or Riemannian 
set-up. The main difficulty dealt with here lies in the lack of geometric structure on the spacetime 
which makes it difficult to deploy the analytical techniques used when a metric is prescribed. 
To circumvent this problem, we assume the existence of a suitable family of mollifiers, which 
are adapted to the sole structure available on the manifold, namely the family of flux fields of 
the conservation law. For this reason, our estimate in Theorem 12.41 is, later on in Theorem 13.21 
specialized to the case that a metric is specified. (See also an earlier result in 1101 .) Importantly, 
all of the estimates established in the present paper are coordinate-free. 

An outline of the paper follows. Section [2] is devoted to the general framework. First of 
all, in Section I2TT1 we recall from [11 J some basic concepts about conservation laws posed on 
a spacetime. We then introduce the notion of admissible mollifiers, which allows us to state 
our main result, in Theorem 12.41 below. Then, in Section l2~4l we discuss our first application 
concerning two conservation laws with distinct flux fields. The rest of the section is devoted to 
the proof of Theorem 12.41 In the second part of the paper encompassing the whole of Section [3). 
we treat the special case of conservation laws posed on a Lorentzian manifold, and state our error 
estimate in Theorem l3.2l We are then in a position to provide two more applications in Section l3~3l 
and !3.4l Section l331 contains a proof of Theorem l3.2l 



2 Error estimates for a spacetime 

2.1 Conservation laws based on differential forms 

Let M be an (N + l)-dimensional manifold (with smooth topological structure), which we refer to 
as a spacetime. Denote by A k (M) the space of all smooth fields of ^-differential forms on M, and 
by d : A ,C (M) — > A k+1 (M) the exterior derivative operator. 

Following LeFloch and Okutmustur [ 1 1 J, we recall the formulation of the initial value problem 
for hyperbolic conservation laws posed on a spacetime. The setting is based on differential forms 
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and the conservation law reads 

d(co(u)) = on M, (2.1) 

where u : M — » ]R is the unknown function and the given family of smooth differential N-forms 
co = co(u) e A N (M) depends smoothly upon the variable u and is referred to as the flux field of the 
conservation law. The field co is said to be geometry compatible if it is exact, that is, 

(dco)(u) = 0, w e R 

Moreover, i2.1) is supplemented with the initial condition 

«k = "o, (2.2) 

where Uq is a given data defined on a hypersurface Jfo C M. Throughout this paper we assume 
that the data Uq and, therefore, the unknown function (thanks to the maximum principle) is 
bounded, that is, 

- C < M < Co, (2.3) 

for some Co > 0. We tacitly restrict all values of u to lie in the interval [-Co, Co], and we point out 
that all of our estimates in this paper depend implicitly on this constant Co- 

Following [11 1, we introduce the following notion of global hyperbolicity: the manifold M is 
foliated by hypersurf aces, 

M={Jx t (2.4) 

f>0 

for some (smooth) time-function t, where each slice !K f has the same topology as 'Kq which we 
assume to be a (smooth), compact, orientable N-manifold. For each T > 0, we use the notation 



M r := [J % t . 



0<t<T 

Moreover, given any hypersurface 'K, we denote by i : J-C — > M the canonical imbedding and by 
i* its pullback, taking forms in A k (M) to forms in A k (H). 

Definition 2.1. Let M be a foliated manifold (|2.4[l and co = co(u) be a flux field. The conservation law 
j2.1\ is said to satisfy the global hyperbolicity condition if for allt > the N-form field i^. a>(0) is a 
volume form on "K t , and there exist constants < c < c independent oft such that for all u e [-Co, Co] 

c_r K d u co{G) < ir^ducoiu) <ci^duCo{G), (2.5) 

as inequalities between N-forms defined on the slice "Kf 

In particular, the condition 02.5b implies that for each non-empty smooth hypersurface e C "Kt 
the integral 

T e d u co(u) 



is positive, and its ratio with J i* e d u co(0) is bounded above and below. 
Recall also the following definition of entropy solutions. 

Definition 2.2. A smooth field of N-forms O = Q(TI) is called a (convex) entropy flux field for the 

conservation law l|2.1|l if there exists a (convex) function U : R — > M such that 

Q(u) = I d u U(v)d u co(v)dv, welR. 
Jo 

A measurable and bounded function u : M — > 1R is called an entropy solution to the Cauchy problem 
I|2.1l l- j2.2|l if for every (compactly supported and smooth) non-negative test-function cp, the following 
entropy inequalities hold for any convex entropy pair (U, O) 

I (dcp AQ(u) + cp(dQ)(u)-cpd u U(u)(dco)(u))+ I (pljc^ Q(u ) > 0. (2.6) 
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Under the global hyperbolicity condition above, the flux of the equation naturally induces an 
((N + l)-dimensional) reference volume form on M, namely 

a := it A d u a>(0). (2.7) 

We emphasize that this volume form depends on, both, the flux field (at the state u = 0) and 
the chosen foliation. On the other hand, a more fundamental structure on M is provided by the 
family of N-form flux co = co(u) which determines the conservation law under consideration. 

Finally, we recall from [11] that, under the global hyperbolicity condition, the initial value 
problem (|2.1| |- ||2.2|I admits a unique entropy solution that depends Lipschitz continuously in the 
L 1 norm upon its initial data. 

2.2 Approximate solutions 

The main objective of this paper is to provide a general framework for the derivation of error 
estimates for hyperbolic conservation laws defined on a spacetime. In our statements and proofs 
we will use the Kruzkov entropies 

U{u,k) = \u-k\, fcelR 

with entropy flux 

Q(u,k) = sgn(w - k)(co(u) - co(k)). 
Let v be an (exact) entropy solution to the conservation law \2.6\ . satisfying therefore 

d(Cl(v, I)) + G(v, I) < 0, kR (2.8) 

in the sense of distributions, where G = G(v, 1) are (N + l)-forms defined by 

G(v, l):=sgn(v-l)(dco)(l). 

Following now Bouchut and Perthame |4] and in order to deal with approximate solutions, it 
is convenient to introduce Radon measures to estimate error terms that arise in an approximate 
version of the entropy inequalities. Furthermore, in our setting, rather than functions or vector 
fields, we have to deal with N-form fields on the manifold M. Hence, we now introduce Radon 
measure-valued fields of N-forms. For instance, a distributional (N + l)-form is an element of the 
dual of the space of test-functions on M, while a (scalar) distribution is an element of the dual of 
the space of (smooth) compactly supported (N + l)-forms. In what follows, to keep the notation 
simple we write, for instance, 

I (PP 

Jm 

for the duality bracket between a Radon-measure (N + l)-form field p and a (continuous) test- 
function (p. 

We are now in a position to write the approximate entropy inequalities satisfied by some 
approximate solution u : M — > M, that is, 

d(Q(u, k)) + G(u, k) < dH k + K k =: E k/ ke~R, (2.9) 

where H k is a family of locally Radon measure-valued N-form fields and K k a family of locally 
Radon measure-valued (N + l)-form fields. We make the following key assumption on these error 
terms: there exist a non-negative Radon measure-valued N-form field an and a non-negative 
Radon measure-valued (N + l)-form field such that for every 1-form test-field y and every 
test-function cp > 



sup 

keK 



\ H k Ay 


< 


| och a y 


1 Jm 




Jm 



sup J cpK k < I cpa K . 
fcER Jm Jm 



(2.10) 
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Next, we need to introduce a suitable generalization of Kruzkov and Kuznetsov's mollifiers. 
This is straightforward on a Riemannian manifold, by using the canonical distance function, but 
in the present formalism, we need the following new notion. 

In what follows, if £ = Q(p, q) = t, VA is a function onMxM, then d p Q and dqC denote its 
differentials with respect to the first and second arguments, respectively. We use the notation 
C p for the function p i-> C,(p, q), for q fixed, and we often specify the integration variable in each 
integral to avoid confusion, by adding a subscript to the volume form under consideration. To a 
sequence of functions C, 6 we associate their supports Ejj := supp^ Recalling that a denotes the 

volume form (|2.7|l , we also write |£| := f £ a for the volume of a set E and, therefore, |E*| = J e6 a q . 
We also use the notation f £ a := |E| _1 J £ a. 

Definition 2.3. Fix a non-negative constant A and a smooth 1-form field f> defined on M. A family 
of non-negative (compactly supported and smooth) test-functions (C 6 )6>o defined on MxM is called an 
(A, /3)-admissible family of mollifiers if the following conditions are satisfied: 

1. Unit mass condition: L Cp j? a q = lfor each p e M. 

2. Sup-norm condition: sup^ C, pi] < \E p \~ l for each p e M. 

3. Differential condition: 



rr \d v ? m hy v „\< \ f { \? v 

JJmxm JmJe» 



•p A y Vi q\ 



for each test (2N + l)-form field y. 
4. Symmetry condition: for each u e [-Co, Co] and with y := d u a)(u) 



I I <Pp,q { d V^M A » A "I + d A>A A yi A a v) 

JJmxm 

<A J } (p Viq a p ha q 

JMJEt 



for every bounded function (p : M X M — > 1R + . 

On a Riemannian manifold the distance function allows one to define such a family of test- 
functions. In general, the test-functions should be defined in each application by taking advantage 
of special properties of the given family of approximate solutions (cf. examples below). 

These assumptions arise from natural requirements on the supports E^, volumes \Et\, and test- 
functions C s . They take the proposed form, due to the lack of metric structure on the manifold. For 
instance, Condition lO is a "smallness" condition on the mollifiers, while Condition |[T) replaces 
the unit integral property of standard mollifiers in Euclidian or Riemannian manifolds. 

Condition <j4jl is a symmetry property, enjoyed in the Euclidian space by C, 6 (p,q) = rj({(p,q)), 
which takes then the much simpler form d p C, 6 = -d q C, 5 , when t(p, q) denotes the Euclidian distance 
function. In a Riemannian or Lorentzian setting, a similar property holds, but one needs to use 
parallel transport to compare d p C, and d q Q (see the proof of Theorem 13.21 below) . Our condition 
above is intended to encompass all situations, when no metric is naturally available on the 
manifold. The inequality in Condition (|4) in Definition ^. 3l is motivated by the following formal 
calculation, in wich M = R and C, coincides with a function /] composed with the distance function 
(: 

d p Cy(p) + d q Qy(q) = d p l\(y(p) - y(q)) 

< \r]'{{{p,q))\ su V \y'\e{p,q). 

Now, /] will usually be a standard mollifier, and thus in this example, |/]'| < <5~ 2 < (6| supprjl) -1 . 
Since i{p,q) < 5, the expression above is bounded by | supp sup K \y'\, and Condition (|4) is a 
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generalization of this formal argument, in an integral form. This is also a condition on the size 
of the support of C 5 , since it somehow generalizes the fact that in Riemannian space, €(p, q) < 5 if 
q is in the support of Cp and if those supports are geodesic balls of diameter 5. This assumption 
is necessary since, without a metric on the manifold, it seems difficult to reproduce the above 
argument in a geometric (i.e., coordinate independent) way. 

Note, however, that if one is interested in some particular problem, it is not hard to express the 
constant A in Condition ((4) as a (possibly coordinate dependent) quantity involving derivatives 
of y. 

As for Condition 10 it amounts to a uniform upper bound on the form fields d p C, by a certain 
1-form field |S. It is analogous to the bound |VC 6 | < C5~ N ~ 2 enjoyed by the standard molifiers in 
the Euclidian and Riemannian cases iHl llOl . 

2.3 Statement of the error estimate 

Our main result in the present paper is now stated. Recall that all values u under consideration 
belong to the interval [-Co, Co]. 

Theorem 2.4 (Error estimate for conservation laws on a spacetime). Consider the conservation law 
j2.1\ with flux field co, posed on a spacetime M satisfying the foliated condition (|2.4| | and the global 
hyperbolicity condition \2.5\ for some c,c. Let Q 5 be an (A, ^-admissible family of mollifiers associated 
with some non-negative constant A and 1-form field jS. Consider two functions u(t),v(t) : 'Kt — > ]R that 
belong to l}($Ki)for each t >0 and are right-continuous in t. Assume moreover that v is an (exact) entropy 
solution to the conservation law i2.8i , and that u satisfies the approximate entropy inequalities j2.9i for 
some Hk, Kk, an, cck satisfying the bounds (|2.10l l and such that t^ ocH belongs to L 1 (JC f )/or all t. Then, the 
following L l -type estimate holds for all 5 > and T > 

f iCl(u,v)< f i l Q(u,v)+R 5 [v] + R 5 [co] + R 5 [a], (2.11) 

where 

R b [v] := sup I i*a p I \v p -v q \B q , 

fe(0,T) J lit JEf, 

R 5 [<v] ■■= I A \da> p (v t] ) A« r dco q (v q ) A a p \, 
Jm t Je* 

R 6 [a]:=\ \ lj6Aa H |+ f |fa H | + f a K , 

Jm t Jk um t Jm t 

and B is an (N + l)-form field in q defined by 

Bq := (2c + TA)a q + T sup \d u dco q (u)\. 

u 

A few remarks about the above theorem are in order. First of all, the terms R 5 [v] and R 6 [<x>] 
are expected to tend to zero with 5. For instance, when a metric is prescribed on the manifold, 
the term R 6 [v] is estimated (see Lemma [3T4l below) like in the classical Euclidian case [4J : R 6 [v] < 
CT6 TV(v(0)), provided v has bounded variation. 

Second, under the regularity assumptions on the flux oj, the term R 6 [<x>] is expected to be of 
order 0(6). However, to establish this property, one needs to control the "size" of the sets E*, but 
this cannot be formulated without a notion of distance on M. In contrast, Theorem l3.2l below will 
specialize to the case of a metric on M and on conservation laws based on vector fields, and we 
will see explicitly that R 6 [w] vanishes with 5. 

Finally, note that the quantity i* Q(w, v) can be seen as a measure of the L 1 -norm of the 
difference between u and v. Indeed, in the Euclidian and Riemannian cases, it reduces to |Wp - 
Pp|dVol(5f,). 
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Before discussing some applications of the above theorem, it is interesting to consider the 
special case where the flux field co is "geometry compatible". 



Corollary 2.5. In addition to the assumptions in Theorem |Z41 assume that the flux field co is geometry 
compatible, in the sense that (dco)(Ti) = Ofor each u e M. Then, the following error estimate 

f ?Q{u,v)< f i*Cl(u,v) + R 6 [v] + R b [a], (2.12) 

holds for all T > 0, with R 6 [v] and R 6 [a] defined as in Theorem YLQ 

Furthermore, ifv is sufficiently smooth so that R 6 [v] — » as 5 — > and if u is also an exact entropy 
solution, then the L l -like distance between u and v 



1 i-» I i*Q(u, v) 



is non-increasing in time. 



The second statement in the above corollary is nothing but the contraction property of the 
semi-group of entropy solutions. 

2.4 Application (I). Comparing two conservation laws 

Theorem 12.41 applies to conservation laws with "modified" flux, and allows us to estimate the 
difference between entropy solutions to two distinct conservation laws. Let co and co be two 
geometry-compatible flux fields, and introduce their corresponding Kruzkov entropy flux field 

Cl(v,k) = sgn(z> - k)(co(v) - co(k)), D.(u,k) = sgn(w - k)(co(u) - co(k)). 

The solutions u, v under consideration satisfy the entropy inequalities 

d(Q(v,k))<0, d(Q(u,k))<0, k e R. (2.13) 

In order to avoid unnecessary technicalities, we may assume that the chain rule applies to ex- 
pressions involving the functions u and v (i.e. bounded functions with bounded variation, for 
instance). 

Theorem 2.6. Let u, v be to entropy solutions satisfying A2.13I I for two flux fields co and co, and assume 
that the conditions in Theorem \2A\ hold for both conservation laws. Then, the following two estimates hold: 

1. Ifv is sufficiently regular so that R 6 [v] — > as 5 — > 0, then 

\ m{u,v)< \ ?Q(u,v) + C \ \d u {co - to) A dv\ (2.14) 

for some uniform constant C > 0. 

2. If, moreover, R s [v] < 5R[v],for some constant R[v] independent of b, then 

f m{u,v)< f i*Q(u,v) + c(r[v] f \pAQ(u,a)\\u(p)\\ 

Jm t J. Wo \ Jm t 

+ f tQ{co,^)\u{p)\ 

with 

Q{co,co) = supMu)-aJ(tt)l/l"l- 

u*0 



(2.15) 
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Proof. In order to apply Theorem l2.41 or more precisely its Corollary |2.51 we need to identify the 
structure of the relevant approximate conservation laws. To this end, we write 

d(Q(u,k)) < d(sgn(u - k)Uco - a)(u) - (cu - =: d(y k (u)). 

To show the estimate J2.14L we set K k := d(y k (u)). Using a weak form of the chain rule, we see 
that 

\Kk\ < \d lt {d) - co) A du\. 

Hence, we arrive at the desired estimate (|2.14t , when S — > and after changing the role of u and 
v. 

Second, to establish l|2.15t , we set H k := yk(u). Given an arbitrary 1-form field and following 
arguments in |4J (say that <u(0) = co(0) = for simplicity), we obtain 

\H k Ay\ < C\a H Ay\ 

with an '■= Q(n>, co) \u\. Therefore, the estimate l |2.15t now follows from Corollary |2.5| by minimiz- 
ing over the parameter <5. This completes the proof of Theorem l2.6l □ 

2.5 Derivation of the error estimate 

Proof of Theorem |2~H Let u satisfy the approximate entropy inequalities l|2.9b , and let v satisfy the 
(exact) entropy inequalities J2.8|l . Let cp = (p PjCj be a smooth, compactly supported function on 
MxM. According to d2.6l l and J2.10I L for each k 6 R and q £ M we have for cp > 



- I d p (p A Q p (u p ,k) + I (pG p (u p ,k) 
Jm Jm 

< J (pE k = - I H k A dp(p + J K k <p 
Jm Jm Jm 

< I \a H A d p (p\ + J a K cp, 
Jm Jm 



>m Jm 

thus by taking k — Vq we find 



J dp(p A Qp(u p , v q ) + J (pG p (Up, v q ) 
Jm Jm 



JM JM 

< J \a H Ad p (p\+ I a K (p- 
Jm Jm 

On the other hand, taking / = u p in d2.8l l gives 



(2.16) 



- | d q (p AQ q (Up,v q ) + I (pG q {v q ,u p ) <Q. (2.17) 
Jm Jm 

Since q is the integration variable, the integrals in l|2.17b may be viewed as real-valued functions 
of p. Therefore, we may integrate this function on the manifold, provided a volume form is used. 
Likewise, we may integrate the inequality J2.16I I in q. Choosing the form a from $2.7) . we obtain 



a q dp<p A Qp(w p , v q ) + a, (pG p (u p , v q ) 
m Jm Jm Jm 



and 



< \ a q \ \a H Ad p <p\ 
Jm Jm 

- I Hp I d q (p A Clq(Up, v q ) + I a p I (pG q (v q ,Up) < 0. 
Jm Jm Jm Jm 
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Summing up the above two inequalities and applying Fubini's theorem, we obtain 

- I I d p (p A Q.p(u p , v q ) Aa q + d q <p A 0. q (u v , v q ) A a p 
JJmxm 

+ I I cp(G p (u p , v q ) Al, + G q (v q , Up) A a p ) (2.18) 
JJmxm 

< J I \a H A d p <p| Aa,+fl X A a 9 cp. 
JJmxm 



Let (C )a be an admissible family of mollifiers as in Definition 12.31 and let Xp be a smooth, 
compactly supported real function on M to be specified later. We choose the test-functions 

<Pp,cj - XpZp /Cj > 

which leads us to 

dp(p = C, b dx + xdpC 6 , d q cp = xd q C 5 , 
and so the inequality | |2.18| | becomes 



-a ■ 

JJmxm 



Cdx A 0.p(u p/ v q ) A a 



| I x(d p Q 6 A Q p (up, v q ) A a q + d q C, 6 A Q q (u p , v q ) A a p ) 
JJmxm 



JMxM 

+ I I xQ 6 (Gp(u p/ v q ) A a q + G q (v q/ Up) A a v ) 
JJmxm 

< I I |a H A d p (A< 6 )l A a ? + a K A a^C 5 , 
JJmxm 

which, with obvious notation, has the form 

h-h+h< h- 

The terms I 2 and I 3 . 

We have 

h~ I I X P (^pC 5 a Qp(« p , » 9 ) Aa, + d,C 6 a Q q (u p , v q ) A a p ) 
JJmxm 



MxM 

= 11 Xp ^C 5 A Q q (u p , V p ) A «p 
JJmxm 

+ || XpdcjC 3 A (0 9 (Mp, 17,) - Qq(Mp, Op)) A a p 
J J MxM 

I I jp dpC 5 A (Qp(Mp, o 9 ) - Qp(«p, Z7 p )) A a, 



+ 

JjMxM 

=: + h,2 + h,3> 
with obvious notations, where we have used 



| I A'p^pC 5 A Op(w p , v p ) A a, 
JJmxm 

XpOp(up,v p ) A d v { I Cqa<|) = 0, 



(2.19) 



by Condition |(T} in Definition |2.3l We now analyze the terms in I2. For the first term, note that 

h,i= I Xp« P ( J d q C, 5 AQ q (u p ,v p j) 
Jm JEi 



J x v a v{ \ Q b dQ q (u p ,v p )). 
Jm JEt 



This integration by parts is possible since Q q (u p ,v p ) depends on p only through u and v, i.e., the 
explicit spatial dependence is on q. 
Next, we have 

h,2 = I I I X P (v q - v p )d q C, b A d v Q q (u p , v*) A a p ds, 
Jo JJmxm 

with = si^ + (1 - s)v„, and a similar expression for J2,3- This gives 

^2,2 + 12,3= J II x(^<? - Vp)(dc,C 6 A d v D. q {u p ,v*) A a p 
Jo JJmxm 

+ d p C 6 A d v Qp(u v , v*) A a 9 Jds, 

which, according to Condition lHJl of Definition 12.31 and the definition of Q, yields, for some 
constant A > 0, 



1MJE% 

Putting these estimates together, we obtain 



h,2 + h,3 < I } Ax p \v q - v p \ ocp A a q . 
JmJe^ 



h<~ I Xp« P ( I C, 5 d0. q (u p ,v p )) 
Jm Je 6 p 

I i A^p|^ ? - v p \a p Aa q . 
JmJe*. 



(2.20) 



+ 

JmJe; 

Now, we estimate the term I3 in | |2.19| |. First, note that 

G p (Up, v q ) A a q + G q (v q , u p ) A a p 
= sgn(M p - v q ){dcv p {v q ) AS, - dco q (u p ) A a p ), 

and thus _ _ 

da> p (v q ) A a q — da> q (u p ) A a p 

= dco p (v q ) Aa q - dco q (v q ) A a p 
+ dco q (v q ) A a p - dco q (u p ) A a p =: Ay + Az. 

From Condition 10, that is, C, 5 < \E p \~ x we find immediately 

I X P I C 5 sgn(u p - v q )Ay < R 6 [a>], 
Jm Je;; 

where R 6 [co] is defined as R 5 [co] in the statement of the theorem, but with Xp/\Ep\ instead of l/|Ep|. 
On the other hand, the term A2 gives 

I X P I C 6 sgn(w p - v q )A 2 = - I Xp« P I C 5 dCl q (u p ,v q ). 
Jm Je"„ Jm Je* 



This leaves us with 

-h< I Xpoip I C, 6 dQ q (u v ,v q ) + R 6 [co]. 
Jm Je'I 
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From the last inequality and | |2.20| |, we obtain 

h~h^ I X P «p I C 6 (dQ q (u p ,v q )-dQ q (ii p ,v p )\ 
Jm Je s p 

+ I X Ax P \v q - v p \a p Aa, + R b [a)]. 

JmJe" p 

The first integral is bounded by 

I X P ttp4 sxvp\d u dco q {u)\\v p -v q \, 

Jm Je* u 



so that we finally find 



h-h 



< J Aa p x P \ (a q + sup\d u dco q (u)\)\v p -v q \ (2.21) 

Jm Je" p u 

+ R 5 [co]. 



The terms Ii and I4. 



Let us now turn to l\, the main term in l|2.19|l . As in |4J, we choose the test-functions x — X e to 
be supported on Uo<t<r+i? $U an d constant within any hypersurface J£f, so that x £ is a function of 
f only. Moreover, we arrange that x e ^ 1/ /C e - 1 on U e < f <jJ£ t/ and that d t x e (t) —> 6f=o - <5f = r as 
e — » 0, where Sf =T is a Dirac mass centered at t. Also, we have d^ 6 = d t x e dt. 
We now find 

h = - II C 6 dx e A Qp(u p , o 9 ) A a, 

J J MxM 



I I C d rfx e A Q p (m p , p p ) A a, 

J J MxM 

I I C'^j" A (Qp(Mp, v p ) - Q p (u p , v q j\ A a q . 

J J MxM 



Consider the last integral, and observe that from the non-degeneracy condition $2.5} we find, for 
any positive real function g supported in Mr, 



I g p \dt Ad u Q\ < 2 sup I I g p i*d u Q){u)dt 
Jm u Jo Jm, 

<2c I gpdt A <9«a>(0) = 2c I g p a p . 
Jm Jm 



JM JM 

Thus, using Condition (0, 



IMxM 
< 



I I C 6 rfx e A (Qp(Mp, v v ) - Op(«p, v q )\ A a, 
JJmxm 

I \dx £ A sup ^QplJ' |up - v q \a q 

Jm v Je;1 

< I 2\dtX £ \cdt A d u Q)p{0)A \v p -v q \a q 
Jm Je* 

= I 2\d t x £ \ca p } \v v - v q \a q . 
Jm Je;1 
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Furthermore, we have from Condition (]} in Definition ^. 31 

- I I Zj'dx e A Q p (u p , v v ) A a q = - I dx e A Q p (u p , v p ). 
JJmxm Jm 

Therefore, recalling l|2.19l l and l|2.21)l . we find 

- J dx e A Q p (u p ,v p ) 
Jm 

< I 2\dtX £ \ca p } \v p -v q \a q 
Jm Je* 

J Aa p x e \ (a q + sup \d u dco q \)\v p - v q \ 
Jm Je* u 

6r 



'M 
+ 



(2.22) 



Let us now take e — » 0. First, we have 



-limsup J dx e hQp(u p ,Vp) 
e^o Jm 

= - limsup I d^it) I i*O p (u p ,v p )dt 
< I i*D. p {u p ,v p ) - I i*Q. p (u pi v p ). 

Similarly, and since \d t x e (t)\ — > 6t=o + Sf=r as e — > 0, we find 

limsup I 2\d t x e \ca p A \v p -v q \a q 
e^o Jm Je* 

I Aa p x e A (a q + sup\d u dco q \)\v p -v q \ 
Jm Je* u 



+ 



< R°[v], 

with R 6 [c>] defined as in the statement of the theorem. Also, it is clear that limsup e ^ R 6 [aj] < 
R 5 [co]. 

Finally, we must deal with the term 

h = I I \a H A d p (x% b )\ Aa q + a K A a q x% b ■ 
JJmxm 

We have, using Condition 10 in Definition 12 . 3 1 and | |2.10|| , 

h< J \dx e /\a H \ I + \ I X £ \Pp Aa HU a q 
Jm Jm Jm Je* 

+ f X £ Uk f 

Jm Jm 

Using Condition (]} and taking the limsup as e — > gives J4 < R' 5 [a]. This completes the proof of 
Theorem|2jU □ 
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3 Error estimates for a Lorentzian manifold 



3.1 Conservation laws based on vector fields 

In this section, we derive error estimates for conservation laws posed on Lorentzian manifolds. 
This is motivated by the fact that Theorem 12.41 is greatly improved if the manifold under con- 
sideration has a metric defined on it. To begin with, one does not need to assume the existence 
of a special family of mollifiers, since these are naturally provided by the metric. Second, one 
can introduce second order error terms which allow us to consider more general approximate 
solutions to conservation laws; namely, we obtain an error estimate for very general nonlinear 
diffusion models; see Theorem l3.5l below. 

Let (M,g) be a time-oriented, (N + 1) -dimensional Lorentzian manifold. Here, g is a metric 
with signature (-, +, . . ., +), and we recall that tangent vectors X 6 T p M at a point p e M can 
be separated into timelike vectors (g(X, X) < 0), null vectors (g(X, X) = 0), and spacelike vectors 
(g(X,X) > 0). The manifold is assumed to be time-oriented, so that we can distinguish between 
past-oriented and future-oriented vectors. The Levi-Civita connection associated to g is denoted 
by V and, for instance, allows us to define the divergence operator div ? . Finally, we denote by 
dV g (or dV g (p), to stress the integration variable) the volume measure determined by the metric 
g 

Following (3), a flux-vector on a manifold is defined as a vector field / = f p (u) depending on a 
real parameter u and the conservation law on (M, g) associated with / reads 

div g (/ p (u)) = 0, u : M -» R. (3.1) 

Moreover, the flux-vector / is said to be geometry compatible if 

div g f p (u) = 0, u e R, p e M, (3.2) 

and to be timelike if its w-derivative is a timelike vector field 

g(d u f p (u), d u f p {u)) < 0, p e M, u e R. (3.3) 

We are interested in the initial-value problem associated with (|3.11 . So, we fix a spacelike 
hypersurface Jfo c M and a measurable and bounded function Uq defined on !Kq. Then, we 
search for a function u = u p e L°° (M) satisfying ll3.ll in the distributional sense and such that the 
trace of u on "Kq coincides with Uq, that is 

"|Mo = "o- (3.4) 

It is natural to require that the vectors d u f p (u), which determine the propagation of waves in 
solutions of l|3.1j) , are timelike and future-oriented. Thus, we will assume throughout that the 
flux-vector in equation l l3.ll is timelike, in the sense of d3.3|) . 

As in the previous section, we assume that the manifold M is globally hyperbolic, which in 
this Lorentzian setting means that there exists a foliation of M by spacelike, compact, oriented 
Riemannian hypersurfaces !K t (f e R): 

teK 

Any hypersurface !Kf is referred to as a Cauchy surface in M, while the family CKt (t e R) is called 
an admissible foliation associated with JC{ . The future of the given hypersurface will be denoted by 

t>o 

Moreover, we denote by n l the future-oriented, normal vector field to each CK t , and by g l the 
induced metric. Finally, along J£f, we denote by X f the normal component of a vector field X, 
thus X* := g(X, n f ). In the following, when there is no risk of confusion, we write F(u) instead of 
F p (u). 
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Definition 3.1. A flux F - F p (u) is called a convex entropy flux associated with the conservation law 
d3.lt if there exists a convex function U : R — » R such that 



F p (u) = f d u ll{u') d u f v {u') du', p e M, u e R. 
Jo 

A measurable and bounded function u = u p is called an entropy solution of the conservation law fl3.1|l - 
H3.2I I if the following entropy inequality 

f g(F{u),V g c\>)dV g + f (div g F)(u)4>dV g 
Jm + Jm + 

+ f go(F(u ), n°) <Px dV ga - f U'(u)(dw g f)(u)cbdV g >0 

J. 'Ho JM+ 

holds for all convex entropy flux F = F p (u) and all smooth functions (p > compactly supported in M+. 
In particular, the requirements in the above definition imply the inequality 

div^ (f(m)) - (div ? F)(u) + U'(u)(div g f)(u) < (3.5) 

in the distributional sense. For well-posedness results for the initial value problem l !3.1t - l|3.4t , see 

mum. 

3.2 Statement of the error estimate 

For convenience, we consider a Riemannian metric ~g associated with the Lorentzian metric g. We 
fix the natural one, that is, in local coordinates where the matrix of the metric g is diagonal, we 
set g n :- -gn and g u :- ga, (i = 2, . . .,N+ 1). For instance, this allows us to consider the distance 
function ( g associated with ~g. In particular, the volume form and divergence operator associated 
to g or ~g are the same. Also, we write B p (r) for the geodesic ball centered at p e M with radius r, 
with respect to the metric ~g. 

Since we shall rely on Kruzkov's family of entropies for the statement of our results as well as 
the proofs, we write the conservation law 03.5b with Kruzkov's entropy flux 

F p {u, k) = sgn(u - k)(f p (u) - f p {k)), k e R. (3.6) 

Thus, we are given a bounded measurable function u satisfying an approximate entropy 
inequality, 

div g ( sgn(u p - k)(f p (u p ) - f p {k))) + sgn(w p - k) div g f p {k) 

< div g H k + K k + div gt (a k V gl L k ). 

Here, the error terms H k , K k , L k are defined as follows: for each k e R, H k is a distributional vector 
field, that is, an element of the space of linear functionals from the space of smooth 1-forms and 
taking values in the space of (scalar) distributions on M. Thus, for each y e A 1 (M), (H k , y) = y(H k ) 
is a distribution on M, which we assume to be a Radon measure. 

The terms K k and L k are Radon measures, and a k are continuously differentiable functions 
defined on M. We suppose that H k , K k , and L k satisfy the following uniform bounds (with respect 
to k), 

\7(H k )\ < a H \y\ \K k \<a K , \L k \ < a L , (3.8) 

for some positive Radon measures an, ock, «l on M. Here, if a is a measure, \a\ denotes its variation 
in the measure-theoretic sense, and y" is the vector obtained from the 1-form y by raising indices 
using the metric. We also assume that the functions a k satisfy, for some a a independent of k, 

l«/tUV \>a k L, < a a . 
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Note that due to the presence of a volume form, a measure on M may be seen indifferently as 
an element of the dual of the space of (N + l)-forms or of the dual of the space of test-functions on 
M. Thus, when we write, for instance, J M a^cpdV-^, this denotes the duality between the (scalar) 
distribution Uk and the (N + l)-form (pdV g , and so no regularity is assumed on ax- 
Now, let v denote the exact solution to the conservation law (|3.51 , i.e. for all I e R 

div, ( sgn(p, - l){f q {v q ) - /,(/))) + sgn(p, - I) dw g f q {l) < 0. (3.9) 

Defining 

G p (u,k) = sgn(w - k) div g f p (k), 

the entropy inequalities read 

divj(Fp(w p , k)) + G v (u v , k) < div g H k + K k + div St (fl /c V & L t ) (3.10) 
for the approximate solution u, and 

dw g (F q (v q ,l)) + G q (v q ,l)<0 ( 3 - n ) 

for the exact solution v. 

Our main result in this section is Theorem |3.2| It gives a precise, quantitative estimate of the 
evolution of 

I P{up,Vp)dV g i. 

Note that this quantity is equal to 

I \u p - v p \dV gt 

J-K t 

whenever the flux function / of the equation is such that f f (u) = u for all u. Therefore our 
estimates have the same form as the usual estimates in [4], where the manifold is flat and the time 
evolution trivial. In the general case, we have F l (u, v) = \f{u) - f(v)\ for every u, »eR. Since f is, 
by assumption, a strictly monotone function of u, this quantity provides an equivalent measure 
of the difference between u and v in the L 1 -norm, which takes into account the geometry of the 
manifold and the structure of the time-evolution of the foliation under consideration. 

Theorem 3.2 (Error estimate for conservation laws on a Lorentzian manifold). Let u be a function 
satisfying the approximate entropy inequalities (13.71 1, d3-8l >, and let v be an exact solution satisfying 
J3.9I >. Suppose also that u,v are right-continuous with values in L 1 (Jf t ) and that, for some T > 0, «h is 
right-continuous from [0, T) with values in L 1 (Jf t ). For ueK and p e M, define the constants 



A := sup Lip M f, Ai := sup sup Lip u (V x /), 

p p XeTpM 

1X1=1 

A 2 := sup Lip p (div ? /), A 3 := sup Lip u (div g /). 

U Tp 



(3.12) 



< 



Then, for all 5 > the estimate 

F t (u p ,Vp)dV g T 

f F\u v , v v ) dV g o + C(E 5 + Er + E 6 H + E b K + E 6 L ) 
holds, where C is a constant (which may depend on N) and 

E b := (TAj + TA 3 + A ) sup f f \v p - v q \ dV g dV g >, 

f£(0,T) JM,Jb p (6) 



(3.13) 
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E 6 f := T sup \3it\5A 2 , E 6 H := f a H dV g , + jf f a H dV g ,dt, 

t€(0,T) JjC UM T " J J Wf 

4 : = f f a KdV g >dt, : = i f f a L a a dV g ,dt. 

JO J3i, " JO Jj{, 

3.3 Application (II). The semi-group contraction property 

Our main result implies a key property of the semi-group of entropy solutions. 

Corollary 3.3. In addition to the assumptions ofTheorem \3.2\ suppose also that the flux f is geometry- 
compatible and that v has bounded total variation on each slice J£ t . Then, for every 5 > 0, the estimate 



J P{u p , v p )dV g T 



< \ P{u p , v p )dV g o + CT5 sup TV(v m ) + C (E 6 H + E b K + E 6 L ) 

holds where C is independent of 5, but depends on f , and the error terms are defined as in Theorem \3J\ In 
particular, ifu is also an exact solution, the function 



1 1 — > 



I F'(u p , Vp) dV g , 



is non-increasing. 



This result is an immediate consequence of Theorem 13.21 and of the following result, which 
provides the link between the term (associated with the regularity of the exact solution v) 
and its total variation. Recall that for each p 6 M, the exponential map exp p : £>o(S) — > B p (5) c M 
provides a diffeomorphism between the ball of radius 5 on the tangent space at p, Bq(6) C T p M, 
and the geodesic ball (according to the Riemannian metric g~) B p (<5) around p. (Here, 5 must be 
small enough.) The exponential map provides a local chart by identifying TpM with R N+1 . Also, 
in what follows, we write dp, dq instead of the volume element on M to keep the exposition 
uncluttered and to stress the integration variable. It is also more convenient to write v(p) instead 

Of Vp. 

Lemma 3.4. Let v 6 BV(JC t ), t > be a solution to the conservation law j3.% . For 5 sufficiently small, 
one has 



n 

J'KiJbm 



v(p) - v(q)\ dqdp 

(5) 



< C6((l + " Lip ^ lb ° (M) ) TV(W + i||(div g/ K|| LW ), 
where j5 = inf U/P d u f'(u,p) > 0. In particular, if the flux f is geometry-compatible, then 

f / \v(p) - v(q)\ dq dp < C6(l + " Up " { (M) ) TV(v mt ). 
JM,JB p (a) P 

Proof. We may assume that the function v is sufficiently smooth since, by a standard density 
argument, the general result then follows for all functions v with bounded variation. Using the 
exponential map, we write 

f f \v(p)-v(q)\dqdp= \ f \v(exp (0))-v(exp (h))\dhdp. 

Jm,Jb„(5) Jm,Jb (6) 
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Now, consider a partition of unity 1 < i < m subordinate to a covering Q, containing a 
neighborhood of Jf f of radius 5, and write If, = J{ ( n Q;. We have 



f > |i>(exp (0))-z;(exp (h))\dhdp 

Jm,Jb (8) 

= E f Hexp (0))-v(exp (h))\dhdp 

i JU, Jb„(6) 

*f E f Wv(exp p (0))-v(exp (h))\dpdh 

Jb (6) ; JUj 

<C(M) sup V I i/>;|i>(exp (0))-z;(exp (h))\dp, 

heB (5); JUi 

where C(M) is a constant depending on the geometry of M. Note that it is not trivial to reverse 
the order of integration above, and it is necessary to use the partition of unity. This is due to the 
fact that there may be no way to globally specify the isomorphisms between B (6) c R N+1 and 
Bo (6) C TM in a smooth way. Indeed, a sufficient condition to be able to do so is that the tangent 
bundle of M is trivial (i.e., diffeomorphic to M X TM) in a neighborhood of 3£ t , which may not be 
the case. Also, 5 must be small enough so that for each p the point exp p (7?) is well defined. We have 
also abused the notation somewhat since "dh" really stands for the determinant of the Jacobian 
of the exponential map, which, by compactness, may be uniformly bounded as a function of p, 
whence the constant C(M). 

Thus, using that exp p (0) = p, we find 

E X ^^ exp p^ ~ ^Pp^))! rf p - E X ^' X '^ exp /^' dsdv 

< <5 J ^ J \pi\Vjv(exp p (sh))\dpds 
\VgV(p)\ dp + A, 



l-K t 

with 

r 1 i r 

^>i\Vjv(exp p (sh)) - Vjj7(exp p (0))| dpds. 
Next, split the gradient of v into its time and spatial components, 



= 6 fU 

Jo i Ju, 



( \Vjv(p)\dp< ( \d t v(p)\dp+ ( \V x -v{p)\dp. 

J3tt *J3{.t K/Jtl 

Now we use the conservation law to estimate the temporal gradient in terms of the spatial 
gradient. Consider any system of local coordinates on the leaf "Kt. A simple computation shows 
that the conservation law reads for smooth solutions 

d u f v (v{p))d t v{v) + <d«f P Wp)),^)h +di VgfvWP)) = °- 
Thus, after integrating on Jf t we find 

f \d t v\dp < - f d u f\d t v\dp 

J-Kt P J-K, 

< l|Lip f f \^v\d V + \ f \dw g fMrWP- 

P Jw) s P Jm, 
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Finally, since for smooth v 

1 \Vp(p)\dp = TV(v\x t ), 



the result will be proved if sup^g A - 0(5 2 ), which is straightforward. This completes the 
proof of Lemma l3.4l □ 



3.4 Application (III). A nonlinear diffusion model 

In this section, we apply our results to a nonlinear diffusion model on a Lorentzian manifold. 
For simplicity, we will consider the geometry-compatible case, in which the divergence of the 
flux vanishes. Also, to shorten the presentation, we assume here that u and v are regular enough 
so that a weak form of the chain rule applies, see Jlj where (in the Euclidian case) the required 
regularity is that u, v e £> Vi oc , which is also our case. Following |4|, let <p : 1R — > ]R be a Lipschitz 
continuous function. We consider the nonlinear diffusion equation on M, 

div g (f p (u)) = A g ,(p(u), 

where A g t denotes the Laplace-Beltrami operator on the leaf 3£ t . One can check that if <p is 
non-decreasing, the approximate entropy inequalities are 

div ? ( sgn(u p - k)(f p (u p ) - f p (k))) < Agi \(p(u) - <p(k)\. (3.14) 

We will obtain, if u and V have bounded total variation, an error estimate in -^/Lip <p, which is the 
usual estimate in s/e when (p{u) = eu, and also a finer estimate in which only v is required to have 
bounded total variation, and u is only required to be bounded in L 1 . 

Theorem 3.5. Let T > and assume that u satisfy the approximate entropy inequalities | |3.14|| with a 
Lipschitz continuous nonlinear viscosity <p. Let v be an exact solution of the conservation law j3.9i . 

1. Suppose that TV(v\^ t ) < V, TV(w|j<; f ) < LI for all < t < T. Then there is a constant C > such 
that 

f P(u p ,v p )dV g T < f P{u v , Vp) dVgO + CT V(Lip </>) VU. (3.15) 

2. Suppose that TV(v\^ t ) < V, and that ^ \u\dV g < U. Then, 

f P(u p ,Vp)dV g T< f P(Up,Vp)dV g o+CT(QU) 1/3 V 2/3 , (3.16) 

with Q = sup !I#0 \4>(u) - (p(0)\/\u\. 
Proof. Apply Corollary |3.3| with as the sole error term. We find (in the sense of measures) 

\H k \<Lip u( p\V g u\=:a H . 
Estimating the term R 6 [v] using Lemma I3~4l gives 

I F t (Up,v p )dV g T 

< I P{u v , Vp) dV g o + CT5 V + C-^r Lip cp U, 
and the estimate (13.151 follows by minimizing with respect to 5. 
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To establish | |3.16| |, write L k = \(j)(u) - (p(k)\ - \(p(0) - (p(k)\, and so \L k \ < Q\n\ =: a L - Now we 
apply Corollary |3.3| and Lemma [3~4l to find 



I P{u p ,V v )dV g r 



< f F'(u p ,v p )dV g o+CT5V + C^QU, 



from which the estimate (|3.16t follows by choosing the optimal value of 5. This completes the 
proof of Theorem l3.5l □ 

3.5 Derivation of the error estimate 

We provide here the proof of Theorem l3.2l Let <p be a smooth, compactly supported function on 
MxM. From (|3J), (j^TTO) and l|3lT)> , we have for all Zc, Z e R 

- I d p (p(F p (u p ,k))dV g (p) + I <p M Gp(u p ,k)dVg(p) 
Jm Jm 

< I (div^ H, c + Kjt + div g /(fl)cV^L, c ))(p M dV g (p) 

JM 

< \ \d<p(H k )\dV g (p)+ \ \K k \(pdV g (p)+ \ \L k \div g ,(a k y gt( p)dV g (p) 
Jm Jm Jm 

< f E(<p)dV g (p), 
Jm 



with 

and similarly 



E(cp) := a H \V g (p\ + a K <p + a L a a (\A g ,(p\ + |V^<p|), (3.17) 



f d q cp(F q (v q ,l))dV g (cj) + f (p P , q G q (v q ,l)dV g (q)<0, 
Jm Jm 



)m Jm 
where we recall that d p (resp., d q ) denotes the differential of a function only with respect to p (resp., 
q). Taking k — v q ,l = u p , integrating over M, summing the above inequalities, and using Fubini's 
theorem, we find 

- I I \d p <p(F p (u p , v q )) + d q <p{F q (u p , v q ))\ dV g (p) dV g (q) 
JJmxm 

+ I I f p JG p (u p ,v q ) + G q (v q ,u p j)dV g (p)dV g (q) (3.18) 
JJmxm 



JJm 



Note also that 



Now take 



E(<p) dV g (p) dV g (q). 

MxM 



G p {u p ,v q ) + G q (v q ,u p ) 
= sgn(«p - v q )(div g f p (v q ) - div g f q (u p )). 

6 



<Pp,q — Xp <? p ,qr 



where the function E," is of the form E, p q - C, 6 ({ g (p, q)), the function % is to be chosen later, and ( g 
denotes the distance function associated with the Riemannian metric g. Also, we take C, 6 to be a 
standard family of mollifiers with respect to the Riemannian metric (for a precise definition, see 
ifTOl ). Here, we only record the properties which we will need, namely 

j/^=l, |V^|<^p |A f a<^. (3.19) 
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In what follows we omit the superscript 5 when there is no risk of confusion. One can check that 



the function Cj satisfies 



dq(g(X) — (K q ,X)j, dp{j(X) - -(K p ,X)j, 



for all vector fields X, where K p is the unit tangent vector at p to the (unique) geodesic connecting 
p and q. Therefore, for all vector fields X, we have 

d q t{X) = C o i g (K q , X)j, d p ^X) = -C° (g(K p/ X) g . (3.20) 

Also, when there is no risk of confusion, we shall simply write C' instead of C' ° %• We find from 
(l3~18l l 

I I ( - E,dxF p {u,v) - x(dp£,F p (u,v) + d q £,F q (u,vj) 
JJmxm 

- xi sgn(« - v){ div g f q (u) - div g f p (v))) dV g (p) dV g (q) 

< ff E(xQdV g (p)dV g (q) r 
JJmxm 

where u = u p and V = V„. We now write the last inequality as 

h < h + h + h, (3.21) 



with 



h= ff -ZdxF p (u,v)dV g (p)dV g (q), 
JJmxm 

h = I I x(d p £, F p (u, v) + d q £ F q (u, vj) dV g (p) dV g (q), 
JJmxm 

h= sgn(w - v){ div s f q (u) - div^ f p (vj) dV g (p) dV g (q), 

JJmxm 

h= ff E(xZ)dV g (p)dV g (q). 
JJmxm 



The term I 2 

Note that the term I2 vanishes in the "homogenous" Euclidian case. Writing h = J M XpF 2 dV g {p), 
we can write 



2= I d q £,(F q {u p ,v p ))dV g (q) + J d q dF q (u p ,v q ) - F q {u p ,v p ))dV g (q) 
Jm Jm 

+ I dpdFp{u p ,v q )-F p (u p ,v p ))dV g {q) ^ 3 ' 22 ^ 
Jm 



)M 

= h,\ + h,2 + h,3, 



since 



I d v E,{F v (u p , v p ))dV g {q) = 0. 
Jm 

To see this, note that F p (u p , v p ) does not depend on q, and thus 

P f ^WVfa) = Ffy( f ZdV{qj\ = Fd p ,(l) = 0. 
Jm Jm 

We have 

I d q £,(F q (up,Vp))dV g (q) = - I ^div^F^Wp^dV^f/), 
Jm Jm 
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since we can integrate by parts in this term, as V and u do not depend on the integration variable 
q. The remaining terms are estimated using the regularity of the flux and the difference \v p - v q \, 
as follows. From Il3.22b we find 



h,2 + h,3 



= I \ ( v q~ v p )(d q E,(d v F q (u p , 37*)) - d p dd v F p (u v , v*)))dV g (q)ds, 
Jo Jm 



(3.23) 



Jo Jm 

with v* = sv q + (1 - s)v„. From J3.6t , ( 13.20b , we have 

d v E, F p (u, v) + d£ F q (u, v) = -C((K P , F p (u, v)) g + (K q/ F q (u, v))j). (3.24) 

Now, consider the parallel transport of vectors along a curve on M as follows. Let y : [0, to] — » 
M be a smooth curve, and V the covariant derivative operator associated to the Riemannian metric 
g. Given < s < t < to, the parallel transport is the operator x S/ t : T } ,( S )M —* Ty{t)M such that, given 
a vector X y ^ 6 Ty^M, then T s ,tX y ( S ) e T y ( f )M is the unique solution of the differential equation 

VyX = 0, T S/ fXy( s )| f=s = Xy( s ). 

For our purposes, it is more important to note that, conversely, one may recover the covariant 
derivative from the notion of parallel transport, using the following relation, 

Vy<p)X = lim ; . (3.25) 

' h^O h 

Also, the parallel transport enjoys the property of preserving scalar products, that is, for all vector 
fields X, Y defined along the curve y, 

(X 7 (s), Y y(s))g r(s) = <T s ,fX y(s) , T s ,,Y }is) )g ym . (3.26) 

Furthermore, if y is a geodesic curve, then its tangent vector is invariant under parallel transport. 
With the notations above, if p = y(0) and q = y(h), then To,^-K p = K q . Using the above properties of 
the parallel transport, we find 

iv v, u iv vt w u/v ?h,oF q (u,v)-F p (u,v) x 
(K p , F p (u, v))g - (K q/ F q (u, v))g = -h(K p , )_ 

= -h{K p ,y Kp ,F p ,(u,v)) I 

with h = £g(j), q) and p* some point on the geodesic from p to q. We write simply |X| for the norm 
of a vector with respect to the reference Riemannian metric to keep the notations simple. From 
623, (|3T23l and (32Zll, wim d v F instead of F, we deduce 

I' 2 < I £ div g F q (u p ,Vp)dVg(q) 
Jm 

+ f f h\v p -v£CK ri (V K d v F) r (u p ,v*))_\dV g (q)ds. 
Jo Jm g 

Now, since \K\ = 1 and h < 5, 

[ f h\v p - vMCK p ; (V Kp ,dvF) f (u pi v*))_\dV g (q) ds 
Jo Jm g 

<Ai f 5\C\\v p -v q \dV g (q), 
Jm 

where Ai is defined in (13.12ft . Thus 

h<~ div s F q (u p , v p )dV g (q) dV g (p) 

JJmxm 

+ Ai ff X t>\C\\v p -v q \dV g (q)dV g (p). 
JJmxm 



(3.28) 
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The term J 3 

We now turn to the term 1$ in ([3.211 1, and write 

div g f p (v q ) - div g f q (Up) = div g f p (v q ) - div g f q {v q ) + div g f q (v q ) - dw g f q {u p ), 

thus 

h = ~ I I sgn(u p - v q )\ div ? f p (v q ) - div g f q (v q )\dV g (q) dV g (p) 

JJmxm 

-II sgn(u p - div ? /,(^) - div ? / 9 (Mp)|dV^(<?) dV g (p). 

JJmxm 

The first term is bounded by 

JJmxm 

while the second term is simply 

I I X-S div^ F q (u p , v q )dV g (q) dV g (p). 
JJmxm 

Thanks to 4328), this leads to 

\h + h\< | I xZ\ div g F q (u p , v q ) - div g F q (u p , v p )\dV g (q) dV g (p) 
JJmxm 

+ ff A lX b\C\\v p -v q \ + bA 2X ^dV g {q)dV g {p). 
JJmxm 

We have 

I I x£.\div g F q (u p ,v q ) - div g F q (u p ,v p )\dV g (q)dV g (p) 
JJmxm 

< I I XZLip u (div g f)\v p -v q \dV g (q)dV g (p), 
JJmxm 

which gives 

\h + h\< ff x(Ai6\C\+A 3 Z)\v p -v q \ + 6A2XZdV g (q)dV g (p). (3.29) 

The term Ji 

We now treat the main term 1\. We take X to be a function which is compactly supported in 
time and constant along the hypersurfaces J£f and, thus, for all tangent vectors Y, we have 
dxQO - dtXY*- First, we have 

h= ff -td tX F p (u p ,v p )dV g {p)dV g (q) 
JJmxm 

- ff td tX (F p {u p , v q ) - P p {u v , v p j) dV g (p) dV g (q). 
JJmxm 



Now, 



ff & t x(F p (u p , v q ) - P p (u p , v v j) dV g (p) dV g (q) 

JJMXM 

< ff UpJ^\d t x\\v p -v q \dV g (p)dV g (q), 
JJmxm 
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where we have used the fact that, since d lt f > 0, we have \ f(v) — f{u)\ = F^u, v). From the last 
inequality and ( 13.281 1, ( 13.291 1, the inequality \3.2\\ becomes 



ff -Zd tX P v {u p ,v v )dV g ( V )dV g {q) 

JJMxM 

< ( I ^(p, q)\v p - v q \ dV g (p) dV g (q) + | \ b X t A 2 dV g (p) dV g (q) + h, 
JJmxm JJmxm 



(3.30) 



with 

Ox = ^(Ai6|C'| + {A3) + \dtx\ZAo (3.31) 

and where the constants A, are given by | |3.12|| . Consider now the first term of (|3.30b . From the 
properties of the mollifiers £, namely j ^ PiCj dV g (q) = 1 for all p E M, we find 

ff -& tX P p (u p ,v p )dVg( P )dV g fa)= f -d t xP p (up,v p )dV g (p). (3.32) 
JJmxm Jm 

We choose the function X = X e e (0/ 1) (which only depends on the t coordinate and the small 
parameter e > 0), to be identically one if f 6 (e, T), supported on the set 

fe(0,T+e) 

so that its derivative, d tX , is supported in 

U 

(6(0,e)U(r,r+e) 

and satisfying x e — » l(M f ) tE( o, T ) as e — > 0. Therefore, since <9f^ approaches 6m - <5;k t , and in view of 
the regularity assumptions on u and v and of (13. 321 , we get 



limsup \ -d t x e F p (u p ,v p )dV g (p) 

7 (3.33) 
I P p {u p ,v p )dV g {p)- I ^(Hp^JdV^). 



Next, consider the second term in d3.30b . It yields the term Ef, depending on the regularity of 
V, as follows. In view of (13.311 1 and since there exists a constant Cn such that 

£A 3 + AiSIC'l < ^(Ai + Ag), sup < C N/ (3.34) 



we find 



and 



lim sup A I I £\d t X e \\v p - v q \ dV g (q) dV g (p) 

e^0 JJMxM 

< C N A max I \v p - v q \ dV g (q) dV g (p) 

t=Q - T J"K t JB p (b) 

limsup ff x e {ZM+Mb\C\)\v v -v q \dV g {q)dV g {p) 

e->0 JJMxM 

<TC N (A 1+ A 3 )sup f / fy-^ldVjfaJdV^). 

fe0,T Jn,JBJ5) 
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Thus, we have 



lim sup J J Oi(p, q)\v p - v q \ dV g {q) dV g (p) 

e^O JJMxM 



(3.35) 



<C N (rA 1 + rA3 + A )sup f f \v p -v q \dV g (q)dV g (p). 

feO,T Jji,JB p (5) 

Now, considering the last term in l|3.30b , we find 

ff 6A 2 xZdV g (q)dV g (p) = 6A 2 [ X dV g (p) < T5 sup \M,\A 2 . (3.36) 
jjmxm Jm te(o,T) 

Finally, consider the error term I4 in (|3.21t , given by | |3.17|| . First, note that 

dptp = Ux + Xd v i, Vgt<p = xVg*£, \'<P = X^g'Z- 

From the properties of the test-functions 4 and X s , B.19|l , from (|3.34b , and using the regularity of 
an, we find 

limsup J J \dx e \£a H dV g (p)dV g (q) < limsup I \d t x £ \a H dV g (p) 

e^O JJMxM e^O JM 

<Yj j a HdV g (p), 

t=0,T 

limsup J I \d v E\x £ a H dV g (q)dV g (p) 

e^O JJMxM 

< lim sup J X e \ \CJa H dV g (q)dV g (p) 

e^O JM J B,, (6) 

r 1 f IWI 

< limsup X e -J 4^rdV g (q)dV g (p) 

e^O JM °JB p (6) O 

< limsup^ J x £ Uh d V g (p) < - f a H dV g (p), 

e^O JM $ JM T 

limsup II E,x e ®K dV g (p) dV g (q) < I aicdVg(p), 

e^O JJMxM JM T 

limsup II aLa a (\Agt(p\ + \Vgi(p\)dV g (p)dVg(q) 

e^O JJMxM 

1 1 r 

<C(- + ^J a L a a dV g (p). 

The estimate (|3.13l l now follows from the inequalities above and from J3.33|l , l !3.35|l , (|3.36t . This 
completes the proof of Theorem l3.2l 
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